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The role played by the heaviest fragment in partitions of multifragmenting hot nuclei 
is emphasized. Its size/charge distribution (mean value, fluctuations and shape) gives 
information on properties of fragmenting nuclei and on the associated phase transition. 

1. Introduction 

Nuclear multifragmentation was predicted long ago - and studied since the early 
80's. The properties of fragments which are issued from the disintegration of hot 
nuclei are expected to reveal and bring information on a phase transition of the 
liquid-gas type. Such a phase transition is theoretically predicted for nuclear mat- 
ter. Nuclear physicists are however dealing with finite systems. Following the con- 
cepts of statistical physics, a new definition of phase transitions for such systems 
was recentl y pr oposed, showing that specific phase transition signatures could be 
expected. ISMS] Different anc { coherent signals of phase transition have indeed been 
evidenced. It is only with the advent of powerful Ait detectors ^ like INDRA^that 
real advances were made. With such an array in particular the heaviest fragment 
of multifragmentation partitions is well identified in charge and its kinetic energy 
well measured by taking into account pulse- hei ght defect in silicon detectors ^ and 
effect of the delta-rays in CsI(Tl) scintillators. This paper emphasizes the im- 
portance of the heaviest fragment properties in relation with hot fragmenting nuclei 
(excitation energy and freeze-out volume) and with the associated phase transition 
of the liquid-gas type (order parameter and generic signal for finite systems). 

2. Heaviest fragment and partitions 

The static properties of fragments emitted by hot nuclei formed in central (quasi- 
fused systems (QF) from 129 Xe+ nat Sn, 25-50 AMeV) and semi-peripheral colli- 
sions (quasi-projectiles (QP) from 197 Au+ 197 Au, 80 and 100 AMeV) have been 
compared in detail on the excitation energy domain 4-10 AMeV. To do that hot 
nuclei showing, to a certain extent, statistical emission features were selected. For 
central collisions (QF events) one selects complete and compact events in velocity 
space (constraint of flow angle > 60°). For peripheral collisions (QP subevents) the 
selection method applied to quasi-projectiles minimizes the contribution of dynam- 
ical emissions by imposing a compacity of fragments in velocity space. Excitation 
energies of the different hot nuclei produced are calculated using the calorimetry 
procedure (see ^ for details). First, it is observed that both the percentage of 
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Fig. 1. Full squares and open circles stand respectively for QF and QP sources; Top: average 
values (left) and standard deviation (right) of the charge of the biggest fragment vs the excitation 
energy per nucleon. Bottom: evolution of the charge asymmetry - with (left) and without (right) 
the biggest fragment - as a function of the excitation energy per nucleon. 



charge bound in fragments and the percentage of light charged particles participat- 
ing to the total charged product multiplicity are the same for QF and QP sources 
with equal excitation energy. Thus such percentages provide a good estimate of the 
excitation energy of hot nuclei which undergo multifragmentation. 

What about the size/charge of the heaviest fragment of partitions, Z\ ? 

In figure [1] (upper part), the evolutions, with the excitation energy, of its mean 
value and of the associated fluctuations are plotted. The average charge of the heav- 
iest fragment, for a given system, first strongly decreases with increasing excitation 
energy and then tends to level off, due to the fixed lowest charge value for frag- 
ments. So the mean value appears as also mainly governed by excitation energy 
and is largely independent of system size and of production modes (see SI for lim- 
itations). This effect was already observed in 03 for two QF sources with charges 
in the ratio 1.5; its occurrence when comparing QF and QP sources would indicate 
that their excitation energy scales do agree, within 10%. The fluctuations, on the 
contrary, exhibit sizeable differences. In the common energy range, the fluctuations 
of Z\ decrease when the excitation energy increases but they are larger for QP 
sources. In this latter case they show a maximum value aroun d 4 AM eV which is 
in good agreement with systematics reported for QP sources in EH3 anc j seems to 
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correspond to the center of the spinodal region as defined by the divergences of the 
microcanonical heat capacity. | 14 | 15 | Differences relative to the fluctuations of Z\ 
for QF and QP sources were also discussed in 031 and a possible explanation was 
related to different freeze-out volumes by comparison with statistical model (SMM) 
calculations. We shall see in the next section that indeed different freeze-out vol- 
umes have also been estimated from simulations. An overview of all information 
related to fragment charge partition is ob tain ed with the generalized charge asym- 
metry variable calculated event by event. 022 To take into account distributions of 
fragment multiplicities which differ for the two sources, the generalized asymmetry 
(Az) is introduced: Az = az/{{Z)-JMf rag — 1). This observable evolves from 1 for 
asymmetric partitions to for equal size fragment partitions. For the one fragment 
events, mainly present for QP sources, we compute the Az observable by taking 
"as second fragment" the first particle in size hierarchy included in calorimetry. 
In the left bottom part of fig. [TJ the mean evolution with excitation energy of 
the generalized asymmetry is shown. Differences are observed which well illustrate 
how different are the repartitions of Zf rag between fragments for QF and QP mul- 
tifragmenting sources. QP partitions are more asymmetric in the entire common 
excitation energy range. To be sure that this observation does not simply reflect 
the peculiar behaviour of the biggest fragment, the generalized asymmetry is re- 
calculated for partitions Mf rag > 1, and noted Az\{Zi}, by removing Z\ from 
partitions (bottom right panel of fig. [p. The difference between the asymmetry 
values for the two source types persists. 



3. Heaviest fragment and freeze-out volume 

Starting from all the available experimental information of selected QF sources 
produced in central 129 Xe+ nat Sn collisions which undergo multifragmentat ion, a 
simulation was performed to reconstruct freeze-out properties event by event. P 3* | 17 | 
The method requires data with a very high degree of completeness, which is crucial 
for a good estimate of Coulomb energy. The parameters of the simulation were 
fixed in a consistent way including experimental partitions, kinetic properties and 
the related calorimetry. The necessity of intr odu cing a limiting temperature related 
to the vanishing of level density for fragments^^in the simulation was confirmed for 
all incident energies. This naturally leads to a limitation of their excitation energy 
around 3.0-3.5 AMeV as observed in. 032 The experimental and simulated velocity 
spectra for fragments of given charges (Z=6, 11, 18 and 27) are compared in fig. [2] 
for the different beam energies; when the statistics are sufficient the agreement 
is quite remarkable. Finally relative velocities between fragment pairs were also 
compared through reduced relative velocity correlation functions | 20 | 21 | 22 |23J ^ gee 
fig. [3|). In the simulation the fragment emission time is by definition equal to zero 
and correlation functions are consequently only sensitive to the spatial arrangement 
of fragments at break-up and the radial collective energy involved (hole at low 
reduced relative velocity), to source sizes/charges and to excitation energy of the 
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Fig. 2. Comparison between the experimental velocity spectra (full points) of fragments of a 
given charge and the simulated ones (histograms). Each row refers to a different fragment charge: 
starting from the top Z = 6, Z = 11, Z = 18, Z = 27. Each column refers to a d iffer ent beam 
energy: starting from the left 32 AMeV, 39 AMeV, 45 AMeV and 50 AMeV. From. H3 



sources (more or less pronounced bump at v re d— 0.02-0. 03c). Again a reasonable 
agreement is obtained between experimental data and simulations, especially at 
39 and 45 AMeV incident energies, which indicates that the retained method and 
parameters are sufficiently relevant to correctly describe freeze-out topologies and 
properties. 

The major properties of the freeze-out configurations thus derived are the fol- 
lowing: an important increase, from ~20% to ~60%, of the percentage of particles 
present at freeze-out between 32 and 45-50 AMeV incident energies accompanied by 
a weak increase of the freeze-out volume which tends to saturate at high excitation 
energy. Finally, to check the overall physical coherence of the developed approach, a 
detailed comparison with a microcanonical statistical model (MMM) was done. The 
degree of agreement, which was found acceptable, confirms the main results and 
gives confidence in using those reconstructed freeze-out events for further studies 
as it is done in. HQ] 

Estimates of freeze-out volumes for QF sources produced in Xe+Sn collisions 
for incident energies between 32 and 50 AMeV evolve from 3.9 to 5.7 V/Vo, where 
Vq would correspond to the volume of the source at normal density. 

To calibrate the freeze-out volumes for other sources, we use the charge of the 
heaviest fragment < z[ > or the fragment multiplicity < M^!^ g > , normalized to 
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Fig. 3. Comparison between the experimental (full points) and simulated (histograms) reduced 
relative velocity correlation functions for all the fragments. The reduced relative velocity between 
two fragments with charges Z{ and Zj (Zi j >4) is defined as v re d=v re i /(Zi + Zj) 1 /' 2 . Each panel 
refers to a different beam energy: 32 AMeV (top left), 39 AMeV (top right), 45 AMeV (bottom 
left) and 50 AMeV (bottom right). From. El 



the size of the source, as representative of the volume or density at break-up. From 
the four points for QF sources and the additional constraint that Z^* = M f rag =l 
at V/V =l, we obtain two relations V/V = fi(z[ N) ) and V/V = f 2 (M$g), torn 
which we calculate the volumes for QF sources at 25 AMeV and for QP sources. 
The results are plotted in fig. 21 with error bars coming from the difference between 
the two estimates using /i and fa; note that error bars for the QP volumes are 
small up to 7 AMeV, and can not be estimated above, due to the fall of < Mj?2 > 

at high energy (see fig. 5 of So only < (z[ N ^ >) can be used over the whole 
excitation energy range considered and the derived function is the following: 
V/V = exp{2A7 - AA7 < {z[ N) >) + 0.86. 

The volumes of QP sources are smaller than those of QF sources (by about 20% 
on the E* range 5-10 AMeV). This supports the explanation discussed previously 
starting from fluctuations of the charge of the heaviest fragment in a partit ion. 

Zi also presents some specific dynamical properties. As shown in | 24 | 23 | f or Qp 
sources, its average kinetic energy is smaller than that of other fragments with 
the same charge. The effect was observed whatever the fragment multiplicity for 
Xe+Sn between 32 and 50 AMeV and for Gd+U at 36 AMeV. The fragment- 
fragment correlation functions are also different when one of the two fragments is 
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Fig. 4. Fragmentation position in the exci- 
tation energy-freeze-out volume plane. The 
four full squares (QF sources) are taken 
from The open rectangle gives the es- 
timated position (with error bar) for QF 
source at 25 AMeV, and t he open circles 
those for QP sources. From. 



Z\. This observation was connected to the event topology at freeze-out, the heavier 
fragments being systematically closer to the centre of mass than the others. 



4. Heaviest fragment and order parameter 

The recently developed theory of universal scaling law s of or der-parameter fluc- 
tuations provides methods to select order parameters. I 2 5 | 26 | j n ^us framework, 
universal A scaling laws of one of the order parameters, to, should be observed: 
(to) a P(to) = (/.((to - (to))/(to) a ) 

where (to) is the mean value of the distribution P(m). A=l/2 corresponds to small 
fluctuations, ~ ( m )i an d thus to an ordered phase. Conversely A=l occurs 
for the largest fluctuations nature provides, ~ (to) 2 , in a disordered ph ase. 
For models of cluster production there are two possible order parameters: the 
fragment multiplicity in a fragmentation process or the size of the largest frag- 
ment in an aggregation process (clusters are built up from smaller constituents). 
The method was applied to central collision samples (symmetric systems with total 
masses ^73-400 at bombarding energies between 25 and 100 AMeV) ! 27 l 28 l The 
total (charged products) or fragment multiplicity fluctuations do not show any evo- 
lution over the whole data set. Conversely the relationship between the mean value 
and the fluctuation of the size of the largest fragment does change as a function of 
the bombarding energy: A ~l/2 at low energy, and A ~1 for higher bombarding 
energies. The form of the Z max distributions also evolves with bombarding energy: 
it is nearly Gaussian in the A=l/2 regime and exhibits for A=l an asymmetric 
form with a near-exponential tail for large values of the scaling variable (see fig. [SJ . 
This distribution is close to that of the modified Gumbel distribution, th e resem- 
blance increasing with the total mass of the system studied and being nearly perfect 
for the Au+Au data. The Gumbel distribution is the equivalent of the Gaussian 
distribution in the case of extreme values: it is obtained for an observable which is 
an extremum of a large number of random, uncorrelated, microscopic variables. 
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Fig. 5. (a) Z max distributions for central Xe+Sn collisions at 25-39 AMeV bombarding energies, 
scaled according to A scaling equation; the dashed curve is a best fit to scaled data using a 
Gaussian distribution, (b) As (a) but for bombarding energies 39-100 AMeV: the dashed curve is 
a best fit to scaled data using the Gumbel distribution. From. 12ZI 



Within the developed theory, this behaviour indicates, for extensive systems, the 
transition from an ordered phase to a disordered phase in the critical region, the 
fragments being produced following some aggregation scenario. However si mul ations 
for finite systems have been performed in the framework of the Ising model which 
show that the distribution of the heaviest fragment approximately obeys the A=l 
scaling regime even at subcritical densities where no continuous transition takes 
place. The observed behaviour was interpreted as a finite size effect that prevents 
the recognition of the order of a transition in a small system. More recently the 
distribution of the heaviest fragment was analyzed within the lattice gas model ^2 
and it was shown that the most important finite size effect comes from conservation 
laws, the distribution of the order parameter being strongly deformed if a constraint 
is applied (mass conservation) to an observable that is closely correlated to the 
order parameter. Moreover the observation of the A=l scaling regime was indeed 
observed in the critical zone but was also confirmed at subcritical densities inside 
the coexistence region. 



5. Heaviest fragment and first order phase transition 

At a first-order phase transition, the distribution of the order parameter in a finite 
system presen ts a characte ristic bimodal behavior in the canonical or grandcanon- 
ical ensemble. I 32 | 33 | 34 [35| -j^g bimodality comes from an anomalous convexity of 
the underlying microcanonical entropy. It physically corresponds to the simul- 
taneous presence of two different classes of physical states for the same value of the 
control parameter, and can survive at the thermod yna mic limit in a large class of 
physical systems subject to long-range interactions. ^ 

Indeed if one considers a finite system in contact with a reservoir (canonical sam- 
pling) , the value of the extensive variable (order parameter) X may fluctuate as the 
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system explores the phase space. The entropy function F(X) is no more addititive 
due to the fact that surfaces are not negligible for finite systems and the resulting 
equilibrium entropy function has a local convexity. The Maxwell construction is no 
longer valid. The associated distribution at equilibrium is P(X)^exp(S(X)-AX) 
where A is the corresponding Lagrange multiplier. The distribution of X acquires 
a bimodal character (see fig. [6]). In the case of nuclear multifragmentation, we have 



shown in the previous section that the size of the heaviest cluster produced in 
each collision event is an order parameter. A difficulty comes however from the ab- 
sence of a true canonical sorting in the data. The statistical ensembles produced by 
selecting for example fused systems are neither canonical no r m icrocanonical and 
should be better described in terms of the Gaussian ensemble, I^Sl which gives a con- 
tinuous interpolation between canonical and microcanonical ensembles. Recently a 
simple weighting of the probabilities associated to each excitation energy bin for 
quasi-projectile events was proposed to allow the comparison with the canonical en- 
semble. SQl That weighting procedure is used to allow a comparison with canonical 
expectations for QP sources produced in Au + Au collisions at incident energies 
from 60 to 100 AMeV. Then, a double saddle-point approximation is applied to 
extract from the measured data equivalent-canonical distributions. 

In this incident energy regime, a p art o f the cross section corresponds to col- 
lisions with dynamical neck formation 221 W e th us need to make sure that the 
observed change in the fragmentation pattern S3 is not trivially due to a change 
in the size of the QP. After a shape analysis in the center of mass frame 231, only 
events with a total forward detected charge larger than 80% of the Au charge were 
considered (quasi-complete subevents). Two different procedures aiming at select- 
ing events with negligible neck contribution were adopted. In the first one 32 (I) by 
eliminating events where the entrance channel dynamics indu ces a forward emis- 
sion, in the quasi-projectile frame, of the heaviest fragment Zi. SSI For isotropically 
decaying QPs, this procedure does not bias the event sample but only reduces the 
statistics. In a second strategy (II) the reduction of the neck contribution is ob- 




ln P(X) 



X 



Fig. 6. Canonical ensemble of finite sys- 
tems. The bimodal equilibrium distribution 
is given by P(X)~exp(S(X)-\X). The fig- 
ure shows the case when the Lagrange mul- 
tiplier A is equ al t o the slope of the common 
tangent (from EE]). 
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Fig. 7. Upper part: measured distribution 
of the charge of the largest fragment nor- 
malized to the charge of the source de- 
tected in Au+Au collisions at three different 
bombarding energies. Lower part: weighted 
distributions obtained considering the same 
statistics for each excitation energy bin. The 
left (right) side shows distributions obtained 
with the data selection method (I) ((II)). 
From. ED 



taincd by keeping only "compact" events by imposing (i) an upper limit on the 
relative velocity among fragments, and (ii) a QP size constant within 10% (see^^ 
for details) . In both cases fission events were removed. S3 

The results obtained with the two different selection methods are given in fig. [3 
To take into account the small variations of the source size, the charge of the heaviest 
fragment Z\ has been normalized to the source size. After the weighting procedure 
(lower part of the figure) , a bimodal behavior of the largest fragment charge clearly 
emerges in both cases. In particular in the case of selection (II), we can see that the 
weight of the low Z\ component, associated to more fragmented configurations and 
higher deposited energy, increases with the bombarding energy before the weighting 
procedure (upper part of the figure). This difference completely disappears when 
data are weighted, showing the validity of the phase-space dominance hypothesis. 

Thos e we ighted experimental distributions can be fitted with an analytic func- 
tion (see SD for more details) . From the obtained parameter values one can esti- 
mate the latent heat of the transition of the hot heavy nuclei studied (Z^70) as 
AE = 8.1(±0.4) stat (+1.2 - 0.9) syst AMeV. Statistical error was derived from ex- 
perimental statistics and systematic errors from the comparison between the two 
different QP selections. 
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